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Abstract 

The dynamics of polymers in a random smooth flow is investi- 
gated in the framework of the Hookean dumbbell model. The ana- 
lytical expression of the time- dependent probability density function 
of polymer elongation is derived explicitly for a Gaussian, rapidly 
changing flow. When polymers are in the coiled state the pdf reaches 
a stationary state characterized by power-law tails both for small and 
large arguments compared to the equilibrium length. The character- 
istic relaxation time is computed as a function of the Weissenberg 
number. In the stretched state the pdf is unstationary and exhibits 
multiscaling. Numerical simulations for the two-dimensional Navier- 
Stokes flow confirm the relevance of theoretical results obtained for 
the 5-correlated model. 

KEY WORDS: Elastic dumbbell model; Coil-stretch transition; Turbu- 
lent transport; Batchelor-Kraichnan statistical ensemble. 

1 Introduction 

A polymer is a long, repeating chain formed through the linkage of many 
identical smaller molecules called monomers. Understanding the behavior 
of a single polymeric chain passively advected by a turbulent flow is the first 
step in the study of polymer transport. At equilibrium a polymer molecule 
coils itself up into a ball of radius Rq- When placed in a non- homogeneous 
flow, the polymer deforms and stretches; if the number of monomers com- 
posing the molecule is large, the end-to-end distance of the polymer, R, 
can become much greater than the equilibrium size i?o- A sharp transition 
from the coiled state to the strongly stretched one is observed when the 
strain rate exceeds a critical value [1,2]. 

Two competing effects determine the deformation of a polymer: the 
stretching by the velocity gradients and the elastic relaxation towards the 
equilibrium shape. Theoretical and numerical results suggest that a wide 
range of scales exists where the elasticity of the polymer can be described 
by Hooke's law [3]. In the absence of flow the configuration of the polymer 
converges exponentially to the equilibrium shape: R w i?o e _t / T , where r is 
the characteristic time of the slowest oscillation mode. In contrast, a non- 
homogeneous velocity field can stretch the polymer and elongate it. The 
characteristic value of velocity gradients is determined by the maximum 
Lyapunov exponent A, which is the average logarithmic growth rate of 
fluid particle separations (see, e.g., ref. [4]). 

In random flows the ratio of the relaxation time to the characteristic 
stretching one is determined by the Weissenberg number Wi = At. When 
Wi < 1 contraction is faster and the extension of the molecule converges 
definitively to the equilibrium size; the polymer, therefore, is in the coiled 
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state. When Wi > 1 the elastic force is weak and the velocity gradients can 
deform considerably the molecule. In this case the dynamics of the poly- 
mer depends strongly on the properties of the advecting flow. Balkovsky 
et al. have shown that in random flows the threshold Wi = 1 marks the 
coil-stretch transition and polymers are highly stretched when Wi exceeds 
one [5,6]. Only very recently has the coil-stretch transition in random 
flows been observed experimentally by Gerashchenko et al. [7]. It should 
be however noted that there can be particular flows where the coil-stretch 
transition never occurs. This generally happens when a strong rotation pre- 
vents the molecules from getting aligned with the stretching direction [8] . 

When polymers are highly stretched the coupling between elastic and 
kinetic degrees of freedom should be taken into account. The feedback on 
the advecting flow is then modeled by an additional term in the stress tensor 
appearing in the Navier-Stokes equations. The most important effect of the 
back reaction on the velocity field is probably the drag reduction (see, e.g., 
ref. [9]). The passive approach, however, is justified when the stress due to 
polymers is much smaller than the viscous stress, which can be estimated 
as v v X, v v being the kinematic viscosity of the fluid [5,6]. 

We focus on passive polymer transport: the velocity field is given and 
reaction effects on the advecting flow are neglected. We study the coil- 
stretch transition when the extension of polymers is much smaller than 
their maximum length R ma x and the correlation time of the velocity field 
is negligible. The polymer is modeled by a Hookean elastic dumbbell and 
the turbulent advecting flow belongs to the Batchelor-Kraichnan statis- 
tical ensemble [10, 11]. The same model was previously investigated by 
Chertkov [12] and Thiffeault [13]. Chertkov derived an approximate equa- 
tion for the right tail of the probability density function (pdf) of the end-to- 
end separation of the dumbbell [12]. Thiffeault obtained asymptotic results 
for the stationary solution of this equation in three dimensions [13]. We 
derive the Fokker-Planck equation for the complete pdf of the end-to-end 
distance of the polymer in general dimension d. We thus obtain the exact 
stationary form and the finite-time behavior of the pdf of the elongation. 
Finally, we compare theoretical results with numerical simulations for the 
two-dimensional Navier-Stokes flow. 

The rest of this paper is organized as follows. Section 2 is devoted to the 
description of the model. In Section 3 we derive the Fokker-Planck equation 
for the probability density function of the extension R. The stationary 
solution and the finite time pdf are discussed in Section 4 and in Section 
5 respectively. In Section 6 we discuss the results of numerical simulations 
of polymer dynamics in the two-dimensional Navier-Stokes flow. Section 7 
is reserved to conclusions. 
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2 Elastic dumbbell model 



A simple description of the complex structure of polymers is provided by 
Rouse's model, which represents a polymer as a chain of N beads connected 
by Hookean springs [14]. The dynamics can be solved by decomposing 
the motion of the molecule into a set of linear normal modes. A typical 
relaxation time is associated with the amplitude of each mode. In many 
physical applications the dynamics of the polymer is dominated by the 
fundamental mode, whose relaxation time is of the order of where ( 
is the friction coefficient of the polymer and k is the spring constant. 

The elastic dumbbell model is a simplification of Rouse's model. Only 
the fundamental oscillation mode is taken into account and the elasticity 
of the polymer is modeled by a single spring connecting the ends of the 
molecule. Two beads of radius I are joined at their centers by a linear zero- 
length spring of stiffness k. For a chain of N identical monomers of length 
L the spring constant is k = 3KT/(N — 1)L 2 , where K is the Boltzmann 
constant and T the temperature [15]. The positions of the beads are defined 
by the vectors X\ and the two beads represent the ends of the polymer 
and their separation R = x-i — X\ determines the elongation of the molecule 
(fig. ^) . The idea of using a dumbbell to simulate a polymer dates back to 
Kuhn [16]. For the sake of completeness, we review the derivation of the 
governing equation for the bead connector (see, e.g., ref. [17]). 

The hydrodynamic drag force that one bead experiences due to its rel- 
ative motion to the fluid is given by Stokes' law: — C[^i — v ( x iit)}, where 
Xi is the velocity of the bead i and v(xi,t) the velocity of the fluid at 
the position of the bead 1 . The friction coefficient is £ = 6irlrj v , where r] v 
denotes the dynamical viscosity of the solvent. 

The size of the beads is assumed to be small enough for their dynamics 
to be influenced by thermal noise. This is modeled by a Gaussian process 
such that [18] 

&(*)) = = S^Sit-t') = 1,2) (1) 

Owing to molecular noise the equilibrium size of the dumbbell is not zero, 
but it can be estimated as Rq = y/KT/k by considering the elastic energy 
of the polymer, E = fci? 2 /2, and its equilibrium value at temperature T, 
E = KT/2. 

Taking into account the elastic force, the hydrodynamic drag, and the 

1 In general we ought to take into account the perturbation of the velocity field due to 
the presence of the other bead and Stokes' law should be modified as follows: — C\iCi — 
v[xi,t) — v'(xi,t)], v' being the perturbation. To a first approximation we neglect this 
correction (referred to as "hydrodynamic interaction") which would cause the largest 
relaxation time to decrease as a result of the cooperative motion among the beads [17]. 
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thermal noise, we obtain the dynamical equations for the beads: 



mil = -k (sci - sea) - C [±i - v(x 1: t)] + ^/^CKT^ 



(2) 



m;r 2 = -fc (as a - a*) - £ [i 2 - u(sc 2 , *)] + ^2QKT$, 2 



In physical applications the elongation of the polymer, however large it is, 
never reaches the viscous scale of the turbulent flow, where dissipation and 
advection balance [19] . Below this scale velocity fluctuations are attenuated 
by the viscosity and the flow is smooth. Thus, during its evolution the 
polymer always moves in a regular velocity field, which at the scale R can 
be approximated by a uniform gradient flow [10]: 



Neglecting inertial effects 2 , we finally derive from (J2J a stochastic ordinary 
differential equation for the end-to-end separation of the polymer [17]: 



where r = (/2k is the relaxation time and £ = (£2 — £i)/\/2 has the same 
properties as £1 and £2 [see eq. JJJ]. Equation J2J governs the dynamics 
of the elongation of a polymer transported by a non-homogeneous flow: 
the velocity gradients stretch the molecule, which reacts elastically and 
attempts to recover the equilibrium spherical shape. Thermal fluctuations 
of the solvent are modeled by the white noise \/2Rq/t 

It should be noted that the infinitely extensible dumbbell is a good 
approximation of a polymer only when the extension of the molecule is 
much smaller than its maximum length R ma x- For very large elongations 
(R < Rmax) the relaxation time becomes a function of R and the resulting 
elastic force is nonlinear [3]. In more realistic dumbbell models, Hooke's 
law is replaced by a force diverging as the size of the polymer approaches 
Rmax (for the description of the dynamics of a FENE dumbbell transported 
by the Batchelor-Kraichnan flow see refs. [13,21]). Our results thus apply 
to the range R <C Rmax ■ 

3 Fokker- Planck equation 

The probability density function of the elongation, P(R,t), obeys the 
Fokker-Planck equation associated with eq. © (see, e.g., refs. [18, 22]) 

2 The masses of the beads can be eliminated in the limit where the time scale of 
Brownian fluctuations of the velocity, m/£, is much smaller than the relaxation time of 
the dumbbell, (/2k. This corresponds to considering over-damped oscillations [20]. 



v(x 2 ,t) = v(xi,t) + (x 2 - Xi) ■ Vv(t) 
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dt dR\ t ) t dRdR y ' 

where summation over repeated indices is understood. In the passive 
regime the turbulent advecting flow has to be specified. Within the context 
of Kraichnan's model, v is a Gaussian random field with zero mean and 
correlation function [II] 

(v l {x,t)v : >(x + r,t')) = 9/ lj {r)5{t-t') 

The velocity field is by definition statistically homogeneous in space and 
stationary in time. The non-realistic property of the Kraichnan ensem- 
ble is the ^-correlation in time, which defines a stochastic process with no 
memory. Nevertheless, the absence of time correlation in the advecting 
flow yields closed equations for the correlation functions of the transported 
field. The study of this model brought important results in turbulent trans- 
port theory; these can be interpreted as the limit of real behaviors as the 
velocity correlation time tends to zero (see ref. [4] for a review). A com- 
plete definition of Kraichnan's model requires writing the covariance tensor 
^ u (r) explicitly. If we impose incompressibility, smoothness, and statisti- 
cal invariance under parity and rotations, Sf 1 ^ (r) takes the form [23] 

3> l3 (r) = D S l] - Di[(d + l)5 l3 r 2 - 2rV] 

where the constant D is the eddy diffusivity of the flow, D\ determines the 
intensity of the fluctuations of the velocity, and d is the physical dimension 
of the flow [d = 2, 3 in practical applications). Consequently, the two-time 
correlation of the velocity gradients takes the form 

(V V (t) VV(i')) = 2£>i [{d + l)5 lk V l - 5 ij S ki - 5 ll 6 jk ] S(t - t') (5) 

It is worth noticing that when v is a Gaussian white noise the advection 
term in eq. J3J defines a multiplicative stochastic process like the ones 
considered in refs. [24-27]. 

By averaging eq. over realizations of the velocity field, we obtain the 
evolution equation for the average probability density function P{R, t) = 
(P(R,t)) v : 



dP 



P- 



1 -^P=^^- (6) 



The term ^J^v % P) v can be computed by Gaussian integration by parts [28- 
30] and by using eq. ©■ 
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To analyze the coil-stretch transition, we can restrict our study to the 
pdf of the norm of R, p(R,t) = J P(R,t)R d ~ 1 dQ, where dQ denotes in- 
tegration over angular variables. All the differential operators appearing 
in eq. © are scalar; therefore it is easily shown that p(R, t) obeys the 
one-dimensional Fokker-Planck equation 

I— ^K.dM + ^MJM m 

where /Ci and IC2 are the drift and diffusion coefficients respectively: 

}Ci(R) = a(l-q)R+ b ( d ~ 1 ' IC 2 (R)=aR 2 + b (8) 

R 

The drift and diffusion coefficients are time-independent due to the sta- 
tionarity of the advecting velocity field. The parameter A = D\d(d~ 1) co- 
incides with the maximum Lyapunov exponent of the Batchelor-Kraichnan 
model, while A = 2X/d is the variance of the asymptotic distribution of 
Lyapunov exponents and is related to the correlation of velocity gradi- 
ents [31-34] (see also ref. [4]). The constant q depends on the Weissenberg 
number Wi = At: q — d(l — Wi)/Wi. While the other parameters a, b, d 
are always positive, the sign of q is positive for Wi < 1 and negative in the 
opposite case. 

Equation (JJJ) should be supplemented by an initial condition p(R, 0) 
and by suitable boundary conditions on the interval (0, 00). To this aim 
the Fokker-Planck equation may be rewritten in the form of a probability 
conservation equation 

dp dJ 



ft + 8R=° (10) 



where the functional 



J[p(R,t)} = )d(R)p(R,t) - A [JC 2 (R)p(R,t)} 

is the probability current [18,22,35]. The boundary conditions may then be 
expressed in terms of J(R,t). The normalization of the pdf for all times t 
implies necessarily: J(0,t) = lim^^oo J(R,t) Vt > 0. Here we consider a 
little stronger conditions, usually imposed to the solution of the Fokker- 
Planck equation on an infinite domain: we search for a solution such that 
the associated probability current vanishes both in zero and at infinity: 

J(0,t)= lim J(R,t) = Vt>0 (11) 

R — kxj 
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These conditions (called reflecting boundary conditions) correspond to re- 
quiring that there is no flow of probability through the boundaries of the 
domain of definition [18,22,35]. 

When i? is set to zero eq. J7J) reduces to the equation obtained by 
Chertkov by replacing the norm of the end-to-end vector with its maxi- 
mum component [12]. This approximate equation gives the right tail of 
the pdf of the elongation. Starting from the equation derived by Chertkov, 
Thiffeault studied the large- R tail of the stationary pdf in three dimen- 
sions [13]. Our equation JJJ governs the time evolution of the complete 
pdf of the elongation in general dimension d. In the following sections 
we shall compute the exact stationary form and the finite-time behavior 
of the pdf of the elongation for a Hookean dumbbell transported by the 
Batchelor-Kraichnan flow. 



4 Stationary solution 

Since the drift and diffusion coefficients do not depend on time, p(R,t) 
might relax to a stationary distribution po(R) independent of time and 
of the initial condition. Under reflecting boundary conditions (|llf) the 
probability current associated with po(R) should vanish identically. The 
stationary pdf may then be sought in the form [18,22,35] 



MR) 



c 



exp 



dy 



(12) 



where the constant C and the lower integration limit i?i are fixed by the 
normalization condition. 

Equation Q has a stationary solution only when q is positive [q is 
defined by eq. I©]- After inserting definitions JHJ and © in eq. Ijl2(l . we 
obtain: 



Wi i? 2x - {q+d)/2 

1 + 



p Q (R) = N R d - 1 i [ l + ^^J (q>0) 
Nq is the normalization constant (see ref. [36] formula 3.252.2) 

> Wi y/ 2 r(«j*) 



N 



d J 



r(l)r(i; 



and r denotes the Eulerian integral of the second kind. As the extension 
of the polymer approaches zero, the stationary pdf vanishes as 



p (R) = 0{R d - 1 ) 



as fl-> 



(13) 
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Small values of R are indeed expected when the stretching term is negligible 
and eq. 10) reduces to an Ornstein-Uhlenbeck process [18]. The stationary 
pdf has a maximum at R ~ Rq and displays an algebraic tail for large R: 

po(R) = 0{R- 1 ~' 1 ) asR^oo (14) 

This property has been predicted by Balkovsky et al. for a generic random 
velocity field [5] and has been confirmed by numerical simulations for two- 
dimensional [37] and three-dimensional [38] flows, and by experimental 
studies in random flows [7]. If q > (and so Wi < 1), the normalization 
integral J„ Po(R) dR is dominated by small R: most of the polymers have a 
linear size of the order of i?o ( n g-El- O n the contrary, for q < (or Wi > 1) 
the normalization integral does not converge. This means that a stationary 
pdf does not exist: most of the molecules are highly stretched and the 
passive approach is no longer appropriate for the Hookean dumbbell model. 
Thus, q = (Wi = 1) is the threshold for the coil-stretch transition [5-7]. 

The characteristic exponent q defined in eq. JHJ) coincides with the expo- 
nent estimated by Balkovsky et al. in the neighborhood of the transition 
(r ps A Indeed, they considered a quadratic approximation of the 
Cramer function [30] associated with large deviations of the Lyapunov ex- 
ponents of the flow, and for a velocity field with the Batchelor-Kraichnan 
statistics, the Cramer function is exactly quadratic for all A [31]. 

For a fixed positive q only a finite number of moments of the stationary 
pdf converge. This is due to the presence of highly stretched polymers. The 
number of diverging moments grows with decreasing q and large elongations 
become more and more probable. The m-th order moment can be computed 
explicitly as 

,oo /Wi\ T (4±u±) r ( q ~ m ) 

(JH-JT *»«>«- ( T ) V(i r } ) <m<9) 

(15) 

This latter expression is in accordance with the approximation deduced by 
Thiffeault for d = 3 [13]. 

5 Finite-time evolution 

The time-dependent solution of the Fokker-Planck equation can be ob- 
tained by separation of variables [18,22,35]. The solution should then be 
sought in the form p(R, t) = p v {R) e~ vt and the study of eq. Q reduces to 
the eigenvalue problem 

%YvP v = -vp v (16) 
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where Jzf FP is the Fokker-Planck differential operator 



^ FP ' = ~or {ICl{Ry) + dm (/C2(i?) ' } 

The general solution of eq. Q can then be written as an expansion in 
terms of the eigenfunctions p v where the coefficients are fixed by the initial 
condition. 

The eigenvalue equation associated with the Fokker-Planck equation Q) 
can be transformed into a Gauss hypergeometric equation and can therefore 
be solved explicitly. The complete solution of the eigenvalue problem 116|) 
is postponed to appendix 1X1 Here we describe only the main results. 

The form of the spectrum of the operator Jzf FP fixes the time behavior 
of the pdf of R. The spectrum of Jzf FP in general consists of a discrete as 
well as a continuous branch: 

(a) for all q, the continuous branch is the set of all real values v > q 2 A/8; 

(b) the form of the discrete branch depends on the value of the expo- 
nent q: 

- for a fixed positive q, the discrete branch is bounded from above 
by q 2 A/8 and consists of a finite number of levels of the form: 
v n = nA(q — 2n), n < q/4. If < q < 4, the only discrete 
eigenvalue is v = 0, which is associated with the stationary pdf; 
with increasing q more levels appear; 

- for negative q, the spectrum of Jz? FP has not a discrete branch. 
What are the implications for polymer dynamics? 

Coiled state In the coiled state (q > 0) the pdf relaxes exponentially 
to a stationary distribution. For q large enough (q > 4) the typical relax- 
ation time is the inverse of the smallest non-zero eigenvalue of the discrete 
spectrum, i.e. Ti = [A(q — 2)] _1 . If < q < 4, the convergence to the 
stationary pdf is governed by the continuous branch and the characteristic 
time is T2 = 8/(q 2 A). It is worth noticing that the relaxation time scale de- 
pends on the parameter q or equivalently on the Weissenberg number. This 
dependence is linear for small Wi and becomes quadratic as the system ap- 
proaches the transition Wi = 1 (such behavior is typical of many nonlinear 
multiplicative stochastic processes arising in statistical physics [24-27]). 

Stretched state We have already noted that in the stretched state (q < 
0) there is no stationary pdf: as the time increases, the large-i? tail of 
p(R, t) rises, polymers become more and more elongated and the passive 
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description is no longer suitable for the Hookean dumbbell model. The 
time evolution of the probability p(R, t) with initial condition S(R — Rq) is 
shown in fig- El We have have computed p(R, t) by means of the expansion 
in terms of the eigenfunctions p v [see eq. (|29|) of appendix 

The long time behavior of the pdf of the elongation may also be charac- 
terized by the flatness factors (i? 2m ) / (R 2 ) m . From eq. Q the moments of 
the elongation satisfy a hierarchy of ordinary differential equations which, 
for even-order moments, can be easily solved by Laplace transformation 
and matrix inversion (see appendix ITj|l . The long time behavior of (R 2m ) is 
exponential with an m-depcndcnt exponent which cannot be deduced from 
the spectrum of the Fokker-Planck operator (see, e.g., ref. [25]): 

(B?™) = 0(e S2m */ T ) (t->oo) (17) 

with S2m = 2m[(l + 2m/o?)Wi— 1]. Equation (|17|l is in accordance with the 
prediction of Chertkov [12]. For a fixed q the large-order flatness factors 
grow exponentially in time: the pdf of the elongation is intermittent and 
deviates increasingly from a Gaussian distribution. 

6 Polymer dynamics in "real" turbulence 

One might wonder whether the conclusions reached for the (5-correlated 
flow are of any relevance to real flows. Therefore, we present here results 
from the numerical integration of polymer dynamics, eq. J3J), with the flow 
obtained from the two-dimensional Navier-Stokes equations: 

VP 

d t v + v ■ Vu = h u v Av - av + f (18) 

P 

The velocity field is driven by the large-scale forcing / which is modeled 
by a random zero-mean, statistically homogeneous and isotropic, solenoidal 
vector field. It counteracts the dissipation due to viscosity v v and friction 
a and allows one to obtain a statistically steady state, characterized by 
a direct enstrophy cascade toward small scales [39]. The ensuing velocity 
field is smooth and random [40]. Numerical integration of eq. <|18[) has 
been performed with a standard pseudo-spectral code fully dealiased with 
second-order Runge-Kutta time stepping, on a doubly periodic domain of 
size L = 2n. 

In fig.0]we show the probability density function for polymer elongation 
in the coiled state. The pdf of R reaches a stationary state characterized 
by a peak around Rq, and power-law tails both for small and large elonga- 
tions, in agreement with the results obtained for the Batchelor-Kraichnan 
case (see eqs. I|13|l and (|14[) . and fig. El- The inspection of the vector field of 
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polymer end-to-end separations R reveals a strong correlation with struc- 
tures of the velocity field, as displayed in fig. |5j Above the threshold, the 
pdfs become unstationary, as shown in fig. similarly to what predicted 
for the short-correlated flow (see fig. |3J| • 

Clearly, the basic aspects of polymer stretching by random flows are 
well captured by the analysis of the <5-correlated case. 

7 Conclusions 

We have studied the behavior of polymer molecules seeded in a smooth 
random flow by means of an elastic dumbbell model. Asymptotic results 
for this model were already obtained in refs. [5, 6, 12, 13]. The study of the 
(5-correlated flow allows a complete analytical treatment and the derivation 
of the exact pdf of polymer elongation. In the coiled state, i.e. below the 
critical Weissenberg number Wi = 1, the pdf is stationary and characterized 
by power laws both for small and for large elongations compared to the 
equilibrium length. Above the transition to the stretched state, the pdf is 
not stationary anymore and exhibits multiscaling. The same features are 
observed also for polymers in a turbulent two-dimensional Navier-Stokcs 
flow: this points to the genericity of the results obtained for the Batchelor- 
Kraichnan flow. 

As a by-product of our analysis, an important and new result is the 
computation of the typical time of relaxation to the stationary regime in 
the coiled state. Two regimes can be identified depending on the value 
Weissenberg number: 

-l 

for < Wi < eJ/(ef + 4) 

for d/(d + 4) < Wi < 1. 

The relaxation time diverges quadratically approaching the critical Weis- 
senberg number Wi = 1. Therefore, experimental measures near the coil- 
stretch transition turn out to be quite delicate because of the long time 
needed to reach the stationary regime. For instance, a large statistics 
should be collected to evaluate the scaling exponent q when the Weis- 
senberg number is near the critical value. 

The Hookean model is inadequate when polymers become highly elon- 
gated. Our analysis can be extended to nonlinear elastic models which take 
into account the finite extensibility of polymers. The stationary pdf of the 
elongation for a FENE dumbbell in a (5-correlated flow has been derived in 
ref. [21]; the computation of the relaxation time-scale is still an on-going 



T/r 




. d(l - Wi) 
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work. Finally, it would be of great interest to extend the results of the 
present work to more realistic models of polymer dynamics, e.g. the Rouse 
model [17], within the framework of short-correlated advecting flow. 



A Finite-time pdf 

The time-dependent solution of the Fokker-Planck equation can be ob- 
tained by separation of variables through the ansatz: p(R, t) = p v (R) e~ vt 
[18,22,35]. The study of eq. 10 thus reduces to the eigenvalue problem 

&wP v = -vpu (19) 
where J2? FP is the Fokker-Planck differential operator 

i?FP ^~4 (/Cl(jR) ' ) + ^ (/C2(i?) ' ) 

Under reflecting boundary conditions Jz? F p is symmetric and negative semi- 
defined with respect to the scalar product 

/•OO 

(f,g)= f(R)g(R)[ P0 (R)}- 1 dR (20) 



Hence, v is real and non-negative. The spectrum of „Sf FP in general consists 
of a discrete as well as a continuous branch. If the level v = belongs to the 
spectrum, the Fokker-Planck equation has a stationary solution coinciding 
with the eigenfunction associated with v = 0. When reflecting boundary 
conditions are imposed 3 the eigenfunctions p v form an orthogonal set with 
respect to the scalar product l|20|) : 

' , . , , i 1 5 vv i discrete spectrum , „ 

Pv (R)p v ,(R)\ Po (R)}- 1 dR={;; P (21) 

\o(y — v) continuous spectrum 

Under the further assumption that the set {p u } is a complete basis, the 
time-dependent pdf can be expressed as 

p(R,t) =J2 A nPu n (R)e- / " t + { dv A{v) Vv {R)e~ vt (22) 

n J 

where the coefficients of the expansion are fixed by the initial condition 
p(R,0): 

p OO 

A n = p(R,0)p Un (R)\ Po (R)]- 1 dR 



3 Jv(0) = \imR—>oo Jv{R) = 0, Ji/ being the probability current associated with p„. 
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and 



A{v)= p(R^)p v {R)[ PQ {R)]- 1 dR 
Jo 



It is worth noticing that the transition probability p(R, t\p, 0) is the partic- 
ular solution corresponding to the initial condition p(R, 0\p, 0) = S(R — p) 
and can be expanded in the form i|22|) as 

p(R,t\p,0) = 

E[Po(p)]"V n (pK„(«)e-^+ j dv[p {p)]- l p v {p)p v (R)e- vt 

The eigenvalue equation associated with the Fokker-Planck equation J7J) 
is 4 



(aR 2 + b) 



d 2 p„ 
dR 2 



a(q + 3)R- 



b(d-l) 
R 



dp v 
dR 



b(d-l) 
a(g + l)+ R2 



P, = (23) 



By the change of variable z — — ei? 2 , e = a/b, this latter equation can 
be transformed into a Gauss hypergeometric equation for the new func- 
tion w u (z) = z^ d '' 2 p v (z). Therefore, the solutions of eq. iffify fulfilling 
reflecting boundary conditions take the form 5 



Pv{R) = N v R F (a u , (3 v ;j; —eR 2 ) 



(24) 



where N v is the normalization constant and F (a„, (3 V ] 7; — ei? 2 ) denotes 
the Gauss hypergeometric function with parameters 



1 



a , 1 - I 2 



Av d q 1 

~ ^=2 + 4 + 4 



— 7 = r (25) 
a 2 



From cqs. (|24|) and (|25() it is clear that only d and g determine the form 
of the eigenfunctions p u ; a and 6 are scale factors for the spectrum of Jzf FP 
and for the variable R. The Gauss hypergeometric function is analytic into 
the complex plane with a cut along the positive real axis from 1 to 00 [42] . 
The eigenfunctions (|24|l are then analytic in their entire domain of defini- 
tion (0, 00). 

4 For notational convenience, in this appendix we use the parameters a, b, d, q defined 
in eq. J5J. 

5 The second solution of the Gauss equation scales as R in a neighborhood of the 
origin for odd d and has a logarithmic singularity at R = for even d [41]. Therefore, it 
does not match reflecting boundary conditions for any d. 
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A.l Discrete branch of the spectrum 

The discrete spectrum of the operator Jz? FP is defined by the first of con- 
ditions l|21l) . From the formula expressing the analytic continuation of the 
hypergeometric series in the neighborhood of infinity (see [36] formulas 
9.132.2 and 9.154-5) it is easily shown that p v is normalizablc only when 
7 — ol v is a non-positive integer. Hence, the discrete levels are 

v n = 2an(q — 2n) n < q/A 

For a fixed positive q the set of discrete eigenvalues is bounded from above 
by aq 2 /4 and consists of a finite number of levels. If < q < 4, the 
only discrete eigenvalue is v = 0, associated with the stationary pdf; with 
increasing q more levels appear (see fig. Ell- If 9 is negative, the spectrum 
of „Sf FP has not a discrete branch and the pdf of the elongation does not 
tend to a stationary limit. 

From Euler's formula the eigenfunctions associated with the discrete 
part of the spectrum are rational functions (see [36] formula 9.131 and [41] 
pp. 46-48): 



p Vn (R) = N n R d - 1 (l + ^ 2 



-(d+q)/2 



(- n )m (n - g/2) m (-a/b) r 

rn—0 



^ {d/2) m m\ 



where (A) m denotes the Pochhammer symbol (A) m — A(A — 1) ... (A + 
m — 1), (A)o = 1 and N n is the normalization constant (see [36] formula 
3.194.3): 



_ 2(f) rf/2 r (^) 

Nri. — 



[r(|)r(f)] 1/2 -„,_,„_,. 



n n 

[EE r (l+™+^(!- n -4 



{-l) m+e (-n) m (-n) e (n - g/2) m (n - q/2) t 
(d/2) m (d/2) t m\a 



-1/2 



For large R the eigenfunctions p Un have a power law decay with an exponent 
depending on the order n 

p Vn (R) = OiR- 1 -^ 211 ) as R -> ex) 

A. 2 Continuous branch of the spectrum 

The continuous spectrum consists of all real values v £ (aq 2 /4, oo). The 
eigenvalues can then be written in the form v = aq 2 /4 + 4au 2 , u being 
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a positive real number, and for the sake of simplicity the orthogonality 
condition l|21l) may be replaced by 



p u {R)p ul {R)[p {R)]- 1 dR = S(u - u') 



(27) 



Thus, the eigenfunctions associated with the continuous branch of the spec- 
trum are 



Pu (R)=N u R d - 1 F^ + ^ 



(28) 



where the normalization coefficient can be computed from the asymptot- 
ical expression of the hypergeometric function (see [36] formula 9.132.2 
and [43]): 



r(f + | + m )r(-f + j U ) 



r (2m) 



For large R the eigenfunctions p u are oscillating functions with amplitude 
decaying as i?~ 1_9 / 2 : 



a ,-(d+q)/4 



The phase S u reads 

5 U = arg 



cos(2wlni? + 5 U ) 



as i? — > oo 



r(2* M ) (f) 1 



r(| + j + »u)r(-f + *«) 



Assuming the completeness of the set of eigenfunctions (|2T)|) and 
we conclude that the transition probability of the elongation of an elas- 
tic dumbbell transported by a velocity field with the Batchelor-Kraichnan 
statistics can be written in the form 

9/4-1 

d U bo(p)]-V(pK(i?)e-^ (u) * (29) 



B Moments of the elongation 

From the Fokker-Planck equation J7J we can derive by direct integration a 
hierarchy of ordinary differential equations for the moments of the elonga- 
tion y m (t) = (R rn )(t): 



y m = am(m - q)y m + bm(d + m — 2)j/ m _ 2 



(30) 
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with initial condition y m (0) = y^. Since we are interested in the long 
time intermittency oip{R,t) 1 we restrict attention to even-order moments 
(m = 2i with integer i). In this case the hierarchy of equations (|3(J|) is 
closed by the normalization condition yo = 1 and the 2i-th order moment 
is the solution of a system of 2i linear differential equations. The solution 
can be obtained by Laplace transformation and direct matrix inversion 6 
(see, e.g., ref. [27] for a description of the method). Define the Laplace 
transform of the 2i-th order moment as 

/>OG 

V2i(s) = / e~ st y 2l (t)dt 
Jo 

Equation (|3C)|I can be rewritten in the form 

V2t = ( s ~ 9n)m% - h 2l y 2 i-2 

where 

gi = ai(i — q) and hi — bi(d + i — 2) 
or in matrix form 

v» = X><. ,:>,l! 

3=0 



with 

A t j = (s - g 2l )8 tj ~ h 2i Sij+i 
By matrix inversion we obtain 



1 TT «2fc 



AZ 1 

i<j 
Equation ff3 If) can then be solved for y2% 



j^h 2] ^{s-g 2k ) 



Inverting the Laplace transform by 

rB-\-ioo 



V2i(t) = 7T- / e st y 2i (s)ds 

2 ™ Jb-ioo 



6 Odd-order moments satisfy the same hierarchy of equations, but the problem has an 
infinite number of dimensions. This case, therefore, is more delicate and would require 
studying the convergence of an infinite series [27] . 



1G 



we finally obtain an explicit expression for the 2i-th order moment as a 
finite sum (see ref. [44] for the inversion of eq. (|32|) ~): 



■=3 + 1 



(33) 



j=0 l=j 



where 



his) 



f{s) 



and f{s) = ]]_(s - g 2k ) 



s - 92t 



k=j 



In the stretched state, the exponent gn is an increasing function of i; there- 
fore even-order moments grow asymptotically like e 92it . In the coiled state, 
all convergent moments converge exponentially to the stationary value H15|) 
like e 92 *; divergent moments behave asymptotically as in the stretched 
state. 
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FIGURE CAPTIONS 



Fig. 1: Elastic dumbbell model. The polymer is modeled by two beads 
connected by a linear spring. Each bead experiences the hydrodynamic 
drag force and thermal noise. 

Fig. 2: Batchclor-Kraichnan model: stationary pdf of polymer elonga- 
tion in three-dimensions; linear scale (left) and power-law tails in logarith- 
mic scale (right). With increasing Weissenberg number, polymers get more 
and more elongated. 

Fig. 3: Batchclor-Kraichnan model: pdf of polymer elongation in the 
stretched state (Wi = 1.5, d — 3) in linear scale on the left and logarithmic 
scale on the right. The pdf is intermittent in time and deviates increasingly 
from a Gaussian distribution. 

Fig. 4: Two-dimensional Navier-Stokes flow: the stationary pdf of 
polymer elongation for Wi = 0.4. The straight lines denote the power laws 
R (left) and ii" 1 " 9 with q = 1.32 (right). 

Fig. 5: Left: snapshot of the vorticity field u> — V x v for two- 
dimensional Navier-Stokes turbulence. Right: Polymer end-to-end sepa- 
ration R at the same time as in the left panel. 

Fig. 6: Two-dimensional Navier-Stokes flow: pdf of polymer elongation 
in the stretched state (Wi = 1.6) for different times t/r = 1, 2, 3, 4. 

Fig. 7: Batchclor-Kraichnan model: eigenfunctions associated with the 
discrete part of the spectrum of the Fokker-Planck operator for q = 13, 
a = b = 1, d = 3. 
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Figure 4: 
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